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ABSTRACT 


Nonlinearity introduced in slow-roll inflation will produce weakly non-Gaussian 
CMB temperature fluctuations. We have simulated non-Gaussian large scale GMB 
maps (including COBE-DMR constraints) introducing an additional quadratic term 
in the gravitational potential. The amount of nonlinearity being controlled by the so 
called nonlinear coupling parameter /„;. An analysis based on the Spherical Mexican 
Hat wavelet was applied to these and to the COBE-DMR maps. Skewness values 
obtained at several scales were combined into a Fisher discriminant. Comparison of the 
Fisher discriminant distributions obtained for different nonlinear coupling parameters 
with the COBE-DMK values, sets a constraint of \fni\ < 1100 at the 68% confidence 
level. This new constraint being tighter than the one previously obtained by using the 
bispectrum by Komatsu et al. 


1 INTRODUCTION 

The COBE-DMK data have been found to be compatible 
with Gaussianity by several methods based on real (Torres 
et al. 1995; Kogut et al 1996; Schmalzing & Gorski 1997; 
Novikov et al. 2000), spherical harmonic (Banday et al. 2000; 
Sandvik & Magueijo 2000; Komatsu et al. 2002) and wavelet 
space (Mukherjee, Hobson & Lasenby 2000; Aghanim, Forni 
& Bouchet 2000; Barreiro et al. 2000; Cayon et al. 2001). 
GMB observations by Boomerang, DASI and MAXIMA-I 
(Netterfield et al. 2001; Pryke et al. 2001; Stompor et al. 
2001) seem to indicate that any possible non-Gaussianity 
present in the data will more likely be produced by non¬ 
standard Inflationary models. 

It is important at the moment to study the character¬ 
istics of the non-Gaussianity introduced by the alternative 
theories to standard Inflation. One of the simplest alterna¬ 
tive scenarios that will generate weakly non-Gaussian CMB 
fluctuations is one based on slow-roll inflation. Investigated 
by Salopek & Bond (1990) and Gangui et al. (1994), this 
model generates non-Gaussianity in matter and radiation 
fluctuation fields through features appearing in the inflaton 
potential. In the CMB context, the model has been studied 
by Komatsu & Spergel (2001), Verde et al. (2001) and Ko¬ 
matsu et al. (2002). The gravitational potential includes now 
a quadratic term whith an amplitude regulated by the so 
called non-linear coupling parameter (see Section 2). Verde 
et al. (2001) conclude in their paper that CMB data will be 
more sensitive to the non-Gaussianity generated by these 
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models than observations of high redshift objects. A semi- 
analytical expression of the bispectrum generated by this 
model is presented in Komatsu & Spergel (2001). The mini¬ 
mum non-linear coupling parameter values that will produce 
detectable non-Gaussianity (based on the bispectrum and 
skewness statitics) in different CMB data sets are estimated. 
These values are however larger than the ones expected from 
slow-roll inflation theories. Komatsu et al. (2002) revised the 
constraint impossed by COBE-DMK data on the non-linear 
coupling parameter based on the normalized bispectrum. 
This constraint (slightly larger than the one estimated by 
Komatsu & Spergel 2001), though very weak to say much 
about slow-roll inflation is in any case interesting in itself. 
Deviations from slow-roll inflation might include larger val¬ 
ues of the non-linear coupling parameter generating non- 
Gaussianities that could be detected by CMB experiments. 

We want to asses in this paper which is the minimum 
value of the non-linear coupling parameter that would be 
constraint by the COBE-DMK data based on a wavelet 
space statistic. Comparison of this constraint with the one 
obtained using the normalized bispectrum is an important 
comparison between two methods that will potentially be 
used in the analysis of future CMB data. At the moment, 
simulations of CMB temperature fluctuations generated by 
the slow-roll inflation model here considered are possible at 
large angular scales. The outline of the paper is the fol¬ 
lowing. The simulations used in this work are described in 
Section 2. The method used to constrain the non-linear cou¬ 
pling parameter is based on the Fisher discriminant built on 


© 2001 RAS 


2 


skewness values corresponding to wavelet coefficient maps 
at several scales. These method is presented in Section 3. 
Results are summarized in Section 4. Section 5 is dedicated 
to discussion and conclusions. 


2 SIMULATIONS 

At large scales, the CMB temperature fluctuations are dom¬ 
inated by the Sachs-Wolfe effect and are therefore pro¬ 
portional to the fluctuations of the gravitational potential. 
Weakly non-Gaussian CMB temperature fluctuations will 
be generated by the introduction of a nonlinear term in the 
gravitational potential (?i(n) 

(f){n) = (t)L{n) -I- fni{(t)L{n)- <(t)^{n) >), 

where (pL refers to the linear part of the gravitational poten¬ 
tial, fni is the nonlinear coupling parameter controlling the 
amount of non-Gaussianity introduced and the brackets in¬ 
dicate a volume average. The linear part of the gravitational 
potential gives rise to CMB temperature fluctuations char¬ 
acterised by a power spectrum that in our case is assumed to 
be a CDM power spectrum (computed with the CMBFAST 
software) with D,b = 0.05, fic = 0.3, flu = 0.65, fh, = 0, Ho = 
65 km/s/Mpc and n = 1 and normalized to Qrms = 18/rfc. 

We generate Gaussian as well as non-Gaussian (as ex¬ 
plained above) HEALPix* (Gorski, Hivon & Wandelt 1999) 
projected CMB maps that include the observational con¬ 
straints affecting the COBE-DMK maps. We apply the 
method to the 4-year 53 and 90 GHz COBE-HMR maps af¬ 
ter adding them using inverse-noise-variance weights. These 
maps have a resolution of Nside = 32 corresponding to a 
pixel size of « 110'. Only the pixels outside the extended 
Galactic cut - described in Banday et al (1997) but ex¬ 
plicitly recomputed for the HEALPix sky maps - are con¬ 
sidered in the computation. Best-fit monopole and dipole 
are removed from the pixels under analysis. To follow these 
COBE-DMR observational constraints, the simulated Gaus¬ 
sian and non-Gaussian maps are convolved with a window 
function characterising the COBE-T>MR one as defined in 
Wright et al. (1994). Noise maps, generated with the same 
observation pattern and rms noise levels as the COBE-DMR 
data, are added to them. The extended Galactic cut is also 
impossed to these simulated maps and a best-fit monopole 
and dipole is as well removed from the pixels outside the cut. 
5000 simulated maps are used in the statistical comparison 
described in the next section. 


3 METHOD 

The aim of this work is to obtain the constraint impossed by 
the COBE-DMIl data on the nonlinear coupling parameter 
f„i. We use an approach based on the information provided 
by wavelets at several scales. There are at the moment only 
two wavelets implemented on the sphere that have been used 
in GMB analyses; the Spherical Mexican Hat wavelet (Gayon 
et al. 2001; Martmez-Gonzalez et al. 2001) and the Spherical 
Haar wavelet (Tenorio et al. 1999; Barreiro et al. 2000). The 
former has proved to be more efficient than the Spherical 

* http://www.eso.org/kgorski/healpix/ 


Haar wavelet in detecting certain kind of non-Gaussianity 
as shown in Martmez-Gonzalez et al. (2001). Moreover, the 
non-Gaussianity introduced in the CMB maps by the non¬ 
linear term added to the gravitational potential, is not ex¬ 
pected to be dominated by any prefered direction. Therefore 
a wavelet with spherical symmetry, as the Spherical Mexi¬ 
can Hat, appears to be more appropriate to analyse these 
data. 

Simulated Gaussian and non-Gaussian maps as well as 
the original COBE-Y)MR maps are convolved with Spher¬ 
ical Mexican Hat wavelets with sizes ranging from 2 (220 
arcmins) to 8 pixels. For these maps we compute skewness 
and kurtosis. The non-Gaussianity introduced by the non¬ 
linear term added to the gravitational potential will be vis¬ 
ible when we compare distributions for these statistics for 
the Gaussian and non-Gaussian cases. Moreover, since the 
nonlinear term is small in amplitude and quadratic, we will 
expect the third order cumulant to be larger than any of the 
higher order cumulants. As an example, the skewness and 
kurtosis distributions of Gaussian and non-Gaussian maps 
at 220 arcmins are compared in Figures 1 and 2. The weakly 
non-Gaussianity introduced in the CMB simulations reflects 
in skewed distributions. As expected, the kurtosis seems not 
to be affected by the nonlinear term. We will therefore only 
use the skewness in our final analysis. 

In order to combine all the information we have at the 
seven scales used in the analyses we use the Fisher discrim¬ 
inant. Being x the 7-dimensional data vector including the 
skewness values at the seven scales, the Fisher discriminant 
t{x) is defined as 

t{x) = (/to - fll)^W~^X. 

Where the o subscript refers to the Gaussian case and the 1 
subscript refers to the non-Gaussian cases. W = Vo+Vi gives 
the sum of the covariance matrixes and refer to the 

mean values obtained for the Gaussian and non-Gaussian 
case respectively. This statistic has been previously used 
in other analysis related to detecting non-Gaussianity in 
the CMB (Barreiro & Hobson 2001; Martmez-Gonzalez et 
al. 2001). We establish the COBE-DMR constraint on f„i 
by comparing the COBE-DMR Fisher discriminant value 
with the Fisher discriminant distributions obtained for non- 
Gaussian simulations. 


4 RESULTS 

Fisher discriminant distributions are presented in Figure 3 
for several non-Gaussian cases characterised by the value of 
the nonlinear coupling coefficient. Comparison of Figures 1 
and 3 shows the expected Fisher discriminant result of in¬ 
creasing the power to distinguish two distributions. Vertical 
bars in Figure 3 are drawn at the COBE-DMR Fisher dis¬ 
criminant values for each non-Gaussian case considered. In 
order to obtain a constraint on we calculate the prob¬ 
ability of observing a Fisher discriminant value larger than 
the COBE-DMR one for each of the non-Gaussian cases. 
These probabilities are presented in Table 1. As one can see 
the COBE-DMR data set a constraint on f„i < 1100 at the 
68% confidence level. 

In all previous calculations we have only considered pos¬ 
itive values of the nonlinear coupling coefficient. To study 
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Figure 1. Distribution of skewness values obtained from simulated Gaussian (solid line) and NonGaussian (solid-dotted line) maps 
after convolution with a Spherical Mexican Hat wavelet of size 220 arcmins. Three different NonGaussian cases are presented with 
= 1100,1300,1500 from top to bottom respectively 


Table 1. Probability of having non-Gaussian values of the Fisher discriminant larger(for the positive f^i values) or smaller (for 
the negative f^i values) than the COBE-DMK values 


fnl 

Prob(%) 

1000 

63.0 

1100 

68.7 

-1100 

98.9 

1200 

70.2 

1300 

76.6 

-1300 

99.4 

1400 

78.6 

1500 

81.2 

2000 

93.4 

2500 

97.9 

3000 

99.3 


© 2001 RAS, MNRAS 000, 1-?? 




4 


wv 220 arcmins 



Kurtosis 



Kurtosis 



Figure 2. Distribution of kurtosis values obtained from simulated Gaussian (solid line) and NonGaussian (solid-dotted line) maps 
after convolution with a Spherical Mexican Hat wavelet of size 220 arcmins. Three different NonGaussian cases are presented with 
= 1100,1300,1500 from top to bottom respectively 


the effect of having a negative fni on the Fisher discriminant 
values, we have simulated 5000 non-Gaussian maps with 
fni = —1100, —1300. The distribution of temperature fluctu¬ 
ations in the map would be expected in this case to be shifted 
towards negative values in the same way as the distribution 
of temperature fluctuations is shifted towards positive values 
for fni > 0 (as can be seen from the skewness distributions 
presented in Figure 1). Therefore, we would expect the dis¬ 
tribution of skewness values in the case of /„; < 0 to be 
towards the left of the distribution for /„; = 0. Same be¬ 
haviour is observed for the Fisher discriminant distributions 
as can be seen in Figure 4. The cases for fni = 1100, 1300 
are also included in that figure for comparison. The power 
to distinguish the Gaussian and non-Gaussian distributions 
is approximately the same for the positive and negative 
values of fni- At the 95% confidence level, the power of 
the Fisher discriminant test is 40.7%, 39.1%, 51.5%, 48.6% 
for fni ~ 1100,-1100,1300,-1300 respectively. The same 
“symmetric” behaviour would therefore be expected for the 
any negative fni value in relation to the positive one. As one 


can see from Figure 4 the locus of the COBE-DMR Fisher 
discriminant value does not change much (the probability 
of obtaining Gaussian values larger than the COBE-DMIl 
ones is approximately the same). The probability of hav¬ 
ing non-Gaussian values smaller than COBE-DMR ones is 
98.9% for the case of fni ~ —1100 and 99.4% for the case of 
fni = —1300. The negative values of fni are therefore better 
constrained than the positive ones. We can conclude that 
the COBE-DMIl data set a constraint on the absolute value 
of the nonlinear coupling parameter of \fni\ < 1100 at the 
68% confidence level. 


In order to account for the possible influence of fore¬ 
grounds we have also analysed the COBE-DMIl data after 
foreground removal (using COBE-DIKBF, data). The results 
obtained are consistent with the ones presented above. 
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Figure 3. Distribution of Fisher discriminant values obtained from simulated Gaussian (solid line) and NonGaussian (solid-dotted line). 
The non-Gaussian simulations are generated with nonlinear coupling coefficients = 1000, 1100,1200, 1300,1400, 1500. Vertical bars 
indicate the locus of the COBE-DMR Fisher discriminant for each case. 


5 DISCUSSION AND CONCLUSIONS 

Weakly non-Gaussian CMB temperature fluctuations can be 
generated by the presence of a nonlinear term in the gravi¬ 
tational potential. The amount of nonlinearity can be con¬ 
trolled by the so called nonlinear coupling parameter fni . At 
large scales, the CMB temperature fluctuations are propor¬ 
tional to the gravitational potential. Simulations of weakly 
non-Gaussian CMB temperature fluctuations maps at these 
large scales are therefore reasonably easy to make. 

The aim of this paper is to set a constraint on the value 
of the nonlinear coupling parameter by using the COBE- 
DMR. data. We perform 5000 simulations of Gaussian and 
several non-Gaussian models characterised by the value of 
fni, impossing the COBE-DMR observational constraints. 
With the aid of the Spherical Mexican Hat wavelet, we are 
able to obtain the skewness of these maps at several scales. 
All this information is combined into the Fisher discrimi¬ 
nant. Distributions of this statistic are used to set the con¬ 
straint on fni by comparing the ones obtained from the non- 
Gaussian simulations to the COBE-DMR values. With this 
method we are able to set \fni\ < 1100 at the 68% con¬ 


fidence level. This is a tighter contraint than the one pre¬ 
viously found by Komatsu et al. (2002). They obtain the 
best-fit value of /„;, by fitting the observed COBE-DMR 
bispectrum to a theoretical non-Gaussian model. They esti¬ 
mate statistical uncertainties of this parameter and place a 
68% confidence limit on /„; as \fni \ < 1500 for the extended 
cut. 

A method based on wavelets in which an efficient com¬ 
bination of information from different scales is used, has 
been shown to be very powerful in order to test the non- 
Gaussianity introduced in CMB maps by a nonlinear term 
in the gravitational potential. This work has been carried 
out at large scales at which the Sachs-Wolfe effect domi¬ 
nates. However tighter constraints are expected to be ob¬ 
tained from the analysis of higher resolution maps. 

REFERENCES 

Aghanim, N., Forni, O. & Bouchet, F.R. 2000, astro-ph/0009463 
Banday, A.J., Gorski, K.M., Bennett, C.L., Hinshaw, G., Kogut, 


© 2001 RAS, MNRAS 000, 1-?? 








6 



Fisher discriminant Fisher discriminant 



Fisher discriminant Fisher discriminant 


Figure 4. Distribution of Fisher discriminant values obtained from simulated Gaussian (solid line) and NonGaussian (solid-dotted line). 
The non-Gaussian simulations are generated with nonlinear coupling coefficients f^i = 1100, —1100,1300, —1300. Vertical bars indicate 
the locus of the COBE-DMH Fisher discriminant for each case. 
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